Mon. Not. R. Astron. Soc. 000, dini (2006) Printed 5 February 2008 (MN WF^ style file v2.2) 



Inhomogeneous surface distribution of chemical elements 
in the eclipsing binary ARAur: A new challenge for our 
understanding of HgMn stars^ 

S. Hubrig^f, J. F. Gonzalez^, I. Savanov^, M. Sch611er\ N. Ageorges^, 
§ ; C.R. Cowley^ and B. Wolff^ 

^ ^European Southern Observatory, Casilla 19001, Santiago, Chile 

, ^ Complejo Astronomico El Leoncito, Casilla 467, 5400 San Juan, Argentina 

■ ^ Astrophysical Institute Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany 

■ Department of Astronomy , University of Michigan, Ann Arbor, MI 48109-1042, USA 

^ ^European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany 

^ [ Accepted 2006 Enero 99. Received 2006 Enero 98 
> ■ 

oo ; 

I . ABSTRACT 

'sj" ' We present the results of a high spectral resolution study of the eclipsing binary 

I ARAur. ARAur is the only known eclipsing binary with a HgMn primary star exactly 

, on the ZAMS and a secondary star still contracting towards the ZAMS. We detect for 

the first time in the spectra of the primary star that for many elements the line profiles 
are variable over the rotation period. The strongest profile variations are found for the 
f~i ■ elements Pt, Hg, Sr, Y, Zr, He and Nd, while the line profiles of O, Na, Mg, Si, Ca, 

Oh' Ti, and Fe show only weak distortions over the rotation period. The slight variability 

Q , of He and Y is also confirmed by the study of high resolution spectra of another 

5^ ' HgMn star: a And. A preliminary modelling of the inhomogeneous distribution has 

I been carried out for Sr and Y. Our analysis shows that these elements are very likely 

. concentrated in a fractured ring along the rotational equator. It may be an essential 

' clue for the explanation of the origin of the chemical anomalies in HgMn stars (which 

. ^ \ are very frequently found in binary and multiple systems) that one large fraction of the 

. ring is missing exactly on the surface area which is permanently facing the secondary, 

' and another small one on the almost opposite side. The results presented about the 

5^ , inhomogeneous distribution of various chemical elements over the stellar surface of the 

primary suggest new directions for investigations to solve the question of the origin of 
abundance anomalies in B-type stars with HgMn peculiarity. 
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1 INTRODUCTION 

The origin of the abundance anomalies observed in late B- 
type stars with HgMn peculiarity is still poorly understood. 
Observationally, these stars are characterized by low rota- 
tional velocities and weak or non-detectable magnetic fields. 
Intrinsic photometric variability has been difficult to detect. 
The most distinctive features of their atmospheres are an 
extreme overabundance of Hg (up to 6 dex) and/or Mn 
(up to 3 dex) and a deficiency of He. As more than 2/3 
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of the H gMn stars are known t o belong to spectroscopic 
binaries iHubrig fc Mathv3ll995^ . the variation of spectral 
lines observed in any HgMn star is usually explained to be 
due to the orbital motion of the companion. The aspect 
of inhomogeneous distribution of some chemical elements 
over the surf ace of HgMn stars has b een for the first time 
discussed bv lHubrig fc Mathvsl 1119951) . From the survey of 
HgMn stars in close SBs it was suggested that some chem- 
ical elements might be inhomogeneously distributed on the 
surface, with, in particular, preferential concentration of Hg 
along the equator. In close SB2 systems where the orbital 
plane has a small inclination to the line of sight, a rather 
large overabundance of Hg was found. By contrast, in stars 
with orbits almost perpendicular to the line of sight, mer- 
cury is not observed at all. The first definitively identified 
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spectrum variability which is not caused by the companion 
ha s recently been reporte d for the binary HrMu star a And 
bv'Wah lgren et al.l (|200J^ and lAdelman et alJ (120021 ). They 
suggested that the spectral variations of the Hg II line at 
A 3984 A discovered in high-dispersion spectra are not due 
to the orbital motion of the companion, but produced by the 
combination of the 2.8 day period of rotation of the primary 
and a nonuniform surface distribution of mercury which is 
concentrated in the equ atorial region, in good c orrespon- 
dence with the results of lHubrig fc Mathv3 lll995ll . 

The ZAMS (Zero- Age Main Sequence) eclipsing binary 
ARAur (HD 34364, B9V+B9.5V) with an orbital period 
of 4.13 days at an age of only 4 x 10^ years belong s to 
the Aur OBI association jNordstrom fc Johanserll994ll and 
presents the best case for a study of evolutionary aspects 
of the chemical peculiarity phenomenon. This is the only 
eclipsing binary with HgMn peculiarity known to date. The 
primary and the secondary eclipses are nearly total since its 
orbital inclin ation is 88.5°, and the radii o f both stars are 
almos t equal llNordstrom fc Johansenlll994f) . IChochol et alJ 
l|l988l) discovered a third body in the system. The existence 
of the as yet unseen third star with a mass of at least 0.51 M0 
has been inferred from a light-tim e effect in the observed 
minim a with a period of 25-27 vears. lNordstrom fc JohansenI 
^^9^ studied the parameters of this multiple system in 
detail through an analysis of the available light and radial 
velocity curves. They concluded that the secondary star is 
still contracting towards the ZAMS, while the primary star 
appears to be exactly on the ZAMS. Their data also con- 
firmed the existence of a third body with an orbital period 
of 24 years and an o rbital eccentricity of 0.17. Remarkably, 
iTakeda et al.l ^1979^ reported that ARAur exhibits varia- 
tions in the Hgll A 3983.8 A and YII A 3982.6 A line pro- 
files. Just before the mid-secondary minimum, the Hgll 
line suddenly became the strongest over all the observed 
phases, showing double-line structure. To explain the vari- 
ations of mercury, Takada suggested that the primary star 
has some inho mogeneities such as a cloud, a spot, or stratifi- 
cation. Later, IZverko et al.l il997ft reached a similar conclu- 
sion about the variable profiles of mercury. However, both 
studies used photographic observations at a low spectral res- 
olution of about 30,000 or less and a S/N of ~50. Certainly, a 
more rigorous study using high quality, high-resolution spec- 
troscopy at a very high S /N ratio should provide convincing 
evidence that the Hgll and YII line profiles are variable. 
Therefore, we decided to obtain high quality UVES spec- 
tra of this system to allow a careful study of the previously 
reported anomalous line profile characteristics of these ele- 
ments. 



2 OBSERVATIONS 

Nine spectra of ARAur have been recorded at ESQ with 
the VLT UV- Visual Echelle Spectrograph UVES at UT2 
in 2005. We used the UVES dichroic standard settings 
covering the spectral range from 3030 A to 10,000 A. The 
slit width was set to 0'.'3 for the red arm, corresponding 
to a resolving power of A/AA f» 1.1 X 10^ For the blue 
arm, we used a slit width of 0'.'4 to achieve a resolving 
power of ~ 0.8 x 10^. The spectra have been reduced by 
the UVES pipeline Data Reduction Software (version 2.5; 



Table 1. The log of the observations of ARAur with UVES. The 
signal-to- noise ratio refers to the spectral region around 3990 A. 
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iBallester et aljf200ff) and using IRAF. The signal-to-noise 
ratios of the resulting UVES spectra range from 210 to 380 
per pixel in the one-dimensional spectrum around 3990 A, 
with the maximum S/N of about 450 around 4500 A. The 
summary of our spectroscopic observations of ARAur is 
given in Table 1. 

The theoretical value for synchrono us rotat ion o f 
ARAur is about 22kms~^ jNordstrom fc Johansen<ll994r) . 
The line widths in the spectra of both components are nearly 
identical. Our estimates using eleven Fe II and Ti II lines in 
the region A A 5000 - 5300 A result in u sin i = 22 ± 1 km s" ^ 
This value is consistent with the assumption of synchronous 
rotation. Therefore, the rotation phases have been calcu- 
lated assuming Prot = Poih = 4.13 days. The ephemeris 
have been taken from the most recent study of ARAur 
by ^Ibavrak ct al. ( 2003) . Unfortunately, since our obser- 
vations have been carried out randomly in service mode, no 
spectrum was obtained at the phases of eclipses. Further- 
more, three spectra have been taken around rotation phase 
0.1 and two spectra around phase 0.8. Hence, we have only 
a coverage over six different rotation phases, which is not 
sufficient to perform a reconstruction of the surface abun- 
dance distribution from the rotationally modulated spectral 
line profiles using the Doppler Imaging technique. On the 
other hand, the observations we have obtained so far allow 
us to detect spectral line profile variations, if present, and to 
carry out modelling of spectral line profiles from test images 
with certain abundance distributions. 



3 SPECTRAL VARIABILITY AND 
MODELLING 

The complete set of abundances in the atmospheres of both 
primary and secondary components o f the ARAur sy stem 
can be found in the paper of Khok hlova et alJ (Il995l) . Of 
special interest is the fact that lines of Y, Zr, Nd, Pt, and 
Hg are only present in the spectra of the primary compo- 
nent. Remarkably, inspecting the behavior of spectral lines 
at different rotation phases we discovered that exactly these 
elements show a distinct spectral variability. The lines of a 
few other elements (O, Na, Mg, Si, Ca, Ti, Fe, Sr, and He) 
are present in the spectra of the primary and the secondary 
component, but they appear variable only in the spectra of 
the primary. No noticeable variability of any element has 
been detected in the spectra of the secondary star. 

To study the behavior of the line profiles and to get a 
better idea about the distribution of various elements on the 
surface of the primary star, we decided to compute for a few 
elements, namely Y, Zr, Sr, and Fe, mean line profiles, where 
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Figure 1. Variations of line profiles phased with the rotation period P = 4.13days; a) YII, b) Zrll, c) PtII A4061.7, d) Hgll A3983.9, 
e) SrII, f) He I A 5015.7, g) Ndlll A 4927.5, and h) Fell. The rotational phase increases from top to bottom - see Fig.ITt,. Error bars in 
the lower right corner indicate the standard error of the line profiles. Individual spectra are shifted in vertical direction. 



a set of unblended lines of similar strength is available. We 
note that mean profiles can not be used for Hg and Pt due 
to different isotopic and hyperfine structures of individual 
Hg and Pt lines. We also checked whether the discovered 
variability of spectral lines is not caused by mutual blending 
of lines originating in the individual binary components. 

The problem of analyzing the component spectra in 
double-lined spectroscopic binaries is a tough one, but, 
fortunately, in the past few years several techniques for 
spectral disentangling have been developed. We applied 
the procedure of d ecomp osition described in detail by 
iGonzalez fc Levatd ll2006t) for each observed phase and 
could fully confirm the detected variability of the lines of 
all elements mentioned above. In Fig.Qwe show the behav- 
ior of the line profiles of some elements over the rotation 
period. The Zrll, Ndlll, PtII, and He I lines appear rather 
weak, but still their variations are well noticeable. Interest- 
ingly, while the behavior of the line profiles of YII, PtII, 
Hgll, SrII, and Ndlll is rather similar over the rotation pe- 
riod, the line profiles of Zrll and He I seem to vary with a 
180° phase shift. In Fig.^ we display the mean line profiles 
of Fe II obtained using six unblended lines. Although we see 
very weak distortions in the individual lines over the rotation 
period, the width of the mean line profile suggests that Fe is 
distributed almost homogeneously over the stellar surface. 
The same kind of slight distortions is also observed in the 
01, Nal, Mgll, Sill, Call, and Till line profiles. Moreover, 
the comparison of the observed line profiles with synthetic 
spectra calculated assuming the atmosphere parameters em- 
ployed bv iNordstrom fc_ Joh anscn (1994) (T^ff = 10950 K, 
logg = 4.33) and assuming Vt = 0.5kms~^) reveals that for 
these elements the overall shape of the line profiles slightly 
deviates from the purely rotationally broadened profiles. In 
the spectrum recorded at rotational phase 0.85, the Hgll 
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Figure 2. Radial velocity curves of ARAur. Small filled and 
open circles are the mean stellar radial velocity of the primary 
and secondary components, respectively. In the lower panel we 
show the radial velocity difference of various elements with respect 
to the mean stellar velocity of the primary. Squares present the 
measurements for the Hg II line A 3984, triangles are the results for 
YII, open circles for SrII, pluses for Call, and filled diamonds for 
Ti II lines. The rms of the residuals of the radial velocity curves is 
about 0.5kms~^. The errors in the radial velocities of individual 
elements are of the order of 1 km s~^ for slightly variable elements 
and about 2kms~^ for strongly variable elements. 
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line at A 3984 appears very broad (Fig-Ei)! even broader 
than the mean Hne profiles of Fe II. As the separation of the 
isotope components of this Hgll line is of the same order 
as the Doppler broadening, such a shape of the line profile 
could be interpreted as a signature of the unusual isotopic 
mixture in that rotation phase. However, this explanation 
has to be proven by a detailed analysis of the Hg isotopic 
composition and application of the Doppler Imaging recon- 
struction method. 

In Fig.|5|we show radial velocity curves of AR Aur. For 
a few elements which are distributed inhomogeneously over 
the stellar surface of the primary, we present in the lower 
panel the radial velocity differences (elements Hg, Y, Sr, 
Ca, and Ti). Since most of the line profiles are asymmetric, 
the center-of-gravity method has been used to measure the 
precise line position. However, the measurements are less 
reliable for the Hgll line A 3984, as discussed above. 

In the following we present the first results of our pre- 
liminary modelling of abundance distributions of SrII and 
YII using the SrII line A 4078 and the YII line A 4900. We 
used the so-c alled direct Dop pler Imaging method recently 
presented by IS avanov fc Strassmeier. {.2006.1 . This method 
allows to compute spectral line profiles from test images 
created by varying the local element abundances (or surface 
temperature inhomogeneities) and their location on the stel- 
lar surface. The model atmosphere has been calculated with 
Kurucz's ATLAS9 program i Kurucz. 1993,1 and the local line 
profiles have be en synthesized with th e program SYNTH 
l|piskuno3[l993). The VALD database llKupka et al.lll999l) 
provided the atomic line data necessary for the spectral syn- 
thesis of the local line profiles. The results of modelling two 
almost opposite equatorial spots with radius R~18° located 
at stellar longitudes Zi~84° and Z2~258°, are presented for 
Sr in Fig.|3on the left. Sr was found overabundant in the at- 
mosphere by 0.5 dex and by about 3.0 dex in spots, compared 
to the solar value. For Y we found an overabundance in the 
atmosphere by 2.0 dex and by about 5.0 dex in the spots. 
However, an employment of a partially fractured equatorial 
ring with a width of ~24° (Fig.|21 right) presents a much bet- 
ter consistency between the calculated line profiles and the 
observed ones. Therefore we conclude that the ring struc- 
ture presents a more realistic distribution of Sr and Y on 
the stellar surface. It may be an essential clue for the under- 
standing of the origin of the chemical anomalies in HgMn 
stars (which are very frequently found in binary and mul- 
tiple systems) that one large fraction of the ring is missing 
exactly on the surface area which is permanently facing the 
secondary, and another small one on the almost opposite 
side. It is intriguing that the r esults of Hg mapping o n the 
surface of a And pre sented by IWahlgren et al.l teOOlll and 
lAdelman et al.l (l2002ll are rather similar, showing a concen- 
tration of Hg in a somewhat non-uniform equatorial band. 
From evaluating the changes of the line profiles over the 
rotation period in the spectra of ARAur, it seems that the 
distribution of Hg, Pt, and Nd shows a behavior comparable 
to Sr and Y. 

Since a few years we have at our disposal a few spectra 
of the HgMn star a And, recorded during five consecutive 
nights in 2001. Similar to the observations of ARAur we 
used the UVES dichroic standard settings covering the spec- 
tral range from 3030 A to 10,000 A, giving a S/N greater 
than 400 in the spectral region around 3990 A and more 
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Figure 3. Sr synthetic spectra (top) and Sr maps (bottom) for 
the two models considered: two spots on the surface (left) and a 
fractured ring (right). 

than 500 around 4500 A. From the inspection of these spec- 
tra we see the variability of Hg previously anno unced by 
IWahkren et alJ f200ll) and lAdelman et al.1 ll2002ll . In addi- 
tion, we find that also He I and YII lines appear slightly 
variable. If these elements are concentrated in a fractured 
equatorial ring, similar to that in ARAur, the significantly 
lower variability of Y and He in a And could be explained by 
a different inclination (74°) of the rotational axis to the line 
of sight. The variability of the line profiles A 5015.6 He I and 
A 4884 Y II in the spectra of a And is presented in Fig.^] 

The primary and the secondary components in the 
ARAur system are of considerable interest as they are very 
young and nearly equal-mass stars. iNordstrom fc JohansenI 
(|l99^ suggested that the third star could be a pre-main se- 
quence star. A direct detection of the third component would 
support this scenario and provide an estimate of the mass 
of this component and the inclination between the two or- 
bital planes. We obtained high-resolution images of ARAur 
in November 2005 with NACO at the VLT. The observa- 
tions were carried out in the Ks band using the camera S13. 
However, no companion has been detected. We note that the 
orbital inclination of the wide orbit is unknown. As we are 
presently almost in the phase of conjunction the companion 
would be within the diffraction limit for inclinations below 
50° . Additionally, the third component could be too faint to 
detect in close proximity to the two more massive stars. 



4 DISCUSSION 

The discovery of an inhomogeneous distribution of various 
elements in the atmospheres of AR Aur and a And chal- 
lenges our understanding of the nature of HgMn stars. Using 
FEROS at the ESO 2.2 m telescope we are currently carry- 
ing out a survey of a large sample of HgMn stars to establish 
the spectral variability due to an inhomogeneous elemental 
distribution in the atmospheres of these stars. The obser- 
vations obtained during the last months led already to the 
discovery of a number of HgMn stars with spots (Gonzalez 
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et al. in preparation), proving that the presence of an in- 
homogeneous distribution on the surface of these stars is 
a rather common characteristics and not just a rare phe- 
nomenon found in the atmospheres of very few HgMn stars. 
The re sults of our study modify the claim of Adclman ct al. 
l|2002() that the variability phenomenon in HgMn stars is 
confined to a narrow range of stellar parameters rather than 
characterizing the entire class of HgMn stars. 

Presently we do not have any explanation for the dis- 
covered distribution of Sr and Y in a fractured ring where 
a large fraction of the ring is missing exactly on the surface 
area which is permanently facing the secondary, and another 
small fraction in the almost opposite phase. We believe that 
factors as the presence of a weak tangled magnetic field, 
tidal distortion, or the reflection effect can play a role in 
the development of anomalies in HgMn stars. Although dif- 
fusion due to gravitational settling and radiative levitation 
have been the most popular explanation for the HgMn star 
abundances in the last decades, the accretion of interstel- 
lar material during the pre-main-sequence phase by HgMn 
binary systems seems to be a more convincing choice to ex- 
plain the surface abundances in the presence of magnetic 
fields. Probably, all these mechanisms have to be taken into 
account in future studies of these stars. 

Typically, inhomogeneous chemical abundance distribu- 
tions are observed only on the surface of magnetic chemically 
peculiar stars with large-scale organized magnetic fields. In 
these stars the abundance distribution of certain elements is 
non-uniform and non-symmetric with respect to the rotation 
axis. The presence of magnetic fi elds in the atmospheres o f 
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HgM n stars has been s t udied by^^^^^^J^^^ 
and iHubrig fc Castellil (1200 ll) . iMathvs fc Hubri 
could demonstrate the presence of quadratic magnetic fields 
in two close doubl e-lined systems with HgM n primary stars, 
74 Aqr and x Lup. lHubrig fc CasteUil (|20^ showed for a few 
HgMn stars evidence for a relative magnetic intensification 
of FeH lines produced by different magnetic desaturations 
induced by different Zeeman-split components. As the rel- 



ative intensification is roughly correlated with the strength 
of the magnetic field, it is a powerful tool for detecting mag- 
netic fields which have a complex structure and are diffi- 
cult to detect by polarization measurements. A magnetic 
field of the order of a few hundred Gauss was rec e ntly de - 
tected in four other HgMn stars by IHubrig et a l.' (''2006b|) . 
The sm all sample of HgMn stars for which iHubrig et all 
i2006d) gathered magnetic field measurements with FORS 1 
in spectropolarimetric mode also included a And, with a de- 
tected variable magnetic field of the order of up to 260 G. 
The structure of the measured field in HgMn stars must be, 
however, sufficiently tangled so that it does not produce a 
strong net observable circular polarization signature. A sce- 
nario how a magnetic field can be build up in binary systems 
has been presented some time ago byiHubrig (1993(). As the 
role that magnetic fields play in the development of HgMn 
chemical peculiarities is unclear yet, systematic searches for 
magnetic fields in a larger sample of these stars should be 
conducted and the structure of those magnetic fields needs 
to be determined. 

Finally, we would like to note that since many stars 
with HgMn peculiarity are found in young multiple systems 
and in young associations, further studies of HgMn stars are 
especially important to understand the evolutionary aspects 
of the chemical peculiarity phenomenon. Our recent study 
of the multiplicity of HgMn stars suggests tha t at least 2/3 
of HgMn stars are formed in multiple systems llHubrig et alJ 
l20,06a|) . It is conceivable that HgMn stars have been synchro- 
nized and that their peculiarities were established as early 
as the PMS phase. Both ARAur and a And appear very 
young, with the secondary in the ARAur system in the fi- 
nal stage of contracting towards the ZAMS. a And was also 
found to be located on the ZAMS iHubrig et al...2005i) . 
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